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ABSTRACT 

The diffuse extragalactic 7-ray background in the MeV region is believed to be 
due to photons from radioactivity produced in supernovae throughout the history of 
galaxies in the universe. In particular, 7-ray line emission from the decay chain 56 Ni — > 
56 Co — > 56 Fe provides the dominant photon source flClayton fe Silk 1969Q . Although 
iron synthesis occurs in all types of supernovae, the contribution to the background 
is dominated by Type la events due to their higher photon escape probabilities. 
Estimates of the star formation history in the universe suggest a rapid increase by a 
factor ~ 10 from the present to a redshift z p ~ 1.5, beyond which it either remains 
constant or decreases slowly. Little is known about the cosmological star formation 
history for redshift exceeding z ~ 5. We integrate the observed star formation history 
to determine the Cosmic Gamma-Ray Background (CGB) from the corresponding 
supernova rate history. In addition to 7-rays from short-lived radioactivity in SNIa and 
SNII/Ib/Ic we also calculate the minor contributions from long-lived radioactivities 
( 26 A1, 44 Ti, 60 Co, and electron-positron pair annihilation). The time-integrated 7-ray 
spectrum of model W10HMM (pinto fe Woosley 1988aj pinto fe Woosley 1988b]) 



was used as a template for Type II supernovae, and for SNIa we employ model W7 



(Nomoto et al. 1984). Although progenitor evolution for Type la supernovae is not 
yet fully understood, various arguments suggest delays of order 1—2 Gy between 
star formation and the production of SNIa's. The effect of this delay on the CGB 
is discussed. We emphasize the value of 7-ray observations of the CGB in the MeV 
range as an independent tool for studies of the cosmic star formation history. If the 
delay between star formation and SNIa activity exceeds 1 Gy substantially, and/or the 
peak of the cosmic star formation rate occurs at a redshift much larger than unity, the 
7-ray production of SNIa would be insufficient to explain the observed CGB and a so 
far undiscovered source population would be implied. Alternatively, the cosmic star 
formation rate would have to be higher (by a factor 2-3) than commonly assumed, 
which is in accord with several upward revisions reported in the recent literature. 



Subject headings: cosmology, gamma rays, supernovae 
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1. Introduction 

The study of extragalactic background radiation in various wavelength bands holds the keys 
to many important astrophysical questions. While the cosmic microwave background is an imprint 
of conditions in the early universe, the high-energy background contains information on galaxy 
evolution, stellar explosions, and processes near supermassive black holes in the cores of active 
nuclei. Here we focus on the energy window between 100 keV and 10 MeV, an observationally 
difficult region for which recent analyses of data from COMPTEL (Kappadath et al. 1996), 
HEAO-1 A4 (Kinzer et al. 1997), and SMM (Watanabe et al. 1997) now provide reasonably 
accurate measurements. Below a few hundred keV the observed Cosmic Gamma-Ray Background 
(CGB) is believed to be due to the superposition of unresolved Seyfert galaxies (e.g., Zdziarski 
1996), while for photon energies above 3 MeV blazars are the dominant source population (e.g., 
Sreekumar, Stecker, &; Kappadath 1997). There is no known galaxy type that can fill the 
gap between the Seyfert galaxies and the blazars. However, the CGB around photon energies 

~ lMeV could be partly, or even completely, due to cumulative 7-ray production in supernovae. 
In particular, 7-ray lines from the reaction chain 56 Ni — > 56 Co — > 56 Fe are abundant enough to 



generate a detectable signal flClayton & Silk 19691; IClayton fc Ward 1975| ; [The et al. 19931) 



Although iron synthesis occurs in all SN types, the supernova 7-ray contribution to the 
CGB is dominated by 7-ray escape from the iron synthesis in Type la supernovae (SNI). This 
can be understood from considerations of their rates (R), yields (M, ejected 56 Ni mass), and 
likelihood of photon escape. The rate-yield product satisfies Ri Mi ~ Rn Mn, because Ri < 0.2Rn, 
Mi ~ O.5M , and Mn ~ 0.1M Q . Here we use the label II to represent core collapse events, lumping 
together Type II and Type Ib/c supernovae. The reason that la's dominate the CGB is thus 
ultimately due to the fact that their photon escape probabilities are much higher than those for 
SNII. To determine the SN-induced CGB it is not necessary to use SN light curves, but one must 
distinguish between radioactive nuclei with lifetimes short enough for 7-ray transport to occur in a 
opaque or semi-opaque expanding atmosphere, for which time-integrated light curves are needed, 
and radioactive nuclei whose 7 rays are produced long after the supernova remnant has become 
transparent. In §2 we introduce the basic formalism using monochromatic lines from long-lived 
isotopes. In §4 we discuss the modifications required to treat 7-ray transport in cases of isotopes 
with short lifetimes. We discuss SNI and SNII separately, and emphasize the need to incorporate 
the delay between star formation activity and its associated SNI activity. We then compare the 
CGB estimates to the observations. 



2. Basic Formalism 

The CGB from supernovae is due to a mixture of 7-rays from various radioactivities. 
Short-lived isotopes release photons that are still affected by the expanding supernova envelope, 
and detailed transport calculations are required to determine the emerging spectrum. This 
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strongly affects photons from the decay of 56 Ni, while long-lived isotopes such as 26 Al and 44 Ti 
release photons into an optically thin medium. We begin the discussion of the formalism by 
first considering long-lived nuclei. The ejection of some mass (M e j) of radioactive material then 
generates, assuming I7 per decay, a total photon number 



N Me 1 = 1 2 x 10 53^1 M _ 4) (1) 

A u 



where A is atomic mass number of the radioactive element, u is the atomic mass unit, and M_4 



M e j / 10 4 M is a normalized eject a mass. 



The production rate of photons is proportional to the supernova rate, which we take to be 
proportional to the star formation rate 

Rsn = ZsnM* (2) 

where M* is the star formation rate in MQy -1 , and £s at is measured in units of Mq 1 . Estimates 
of the present-day star formation rate in the Milky Way vary significantly, but within a factor 
two one finds M* ~ few MQy -1 and Rstv ~ few 10~ 2 y _1 (e.g., Tunnies et al. 1997). Therefore, 



we expect that £sn is of order 10 Mq . We assume that this normalization can be universally 
applied to star formation, regardless of galaxy type and whether or not a given galaxy undergoes 
quiescent or bursting star formation. We thus simply apply the scaling relation to the global star 
formation rate density. Recent estimates of this density in the local universe using H Q observations 
( Cole et al. 1994 ; pallego et al. 1995|) yield the global star formation rate density, 



p,{z = 0) = p*(0) ~ 3.7 x lO- 2 h 3 M Mpc~ 3 y" 1 , (3) 

where the Hubble constant is Ho = 100 h km s _1 Mpc -1 . This value is rather similar to earlier 
estimates based on colors of field galaxies ( [Tinsley Danly 1980 ). The corresponding local 7-ray 



emissivity (measured in 7Mpc 3 y x ) is given by 

n 7 (0) = p*(0) &iV N T (4) 

Here we ignore the fact that on scales less than 300 Mpc (redshift less than z = 0.1) the matter 
distribution, and thus also the star formation distribution, is non-uniform. We comment on this 
point briefly in the Conclusions. The 7-ray production is treated as a uniform source density, 
which evolves with redshift in accordance with the evolving star formation rate. 

The emissivity is then integrated over distance (redshift) to yield the differential CGB flux 

— ( 1 cm- 2 s- 1 sr- l ) = F C, (5) 
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where C is a dimensionless combination of cosmological factors, and with the usual definition of 
Hubble length, Lh = c / Ho, the normalization is given as 

F = (47r)- 1 n 7 (0) L H ~ 3.53 x \QT z h 2 A- 1 M^ SN . (6) 

The main integral to be performed is 

POO 

C= / dzil + z^Eiz)- 1 ^) , (7) 

J 

where rj(z) contains evolutionary effects to be discussed below, the factor (l+z) _1 accounts for the 
dilation of the supernova rate, and the function E(z) (equation (13.3) in Peebles 199j ) 



e(z)= n(i + + n R (i + zf + q a 1 1/2 , (8) 



represents the evolution of the Hubble "constant" 



a 



- = H E(z) = H(z) , (9) 
a 

with a(t) being the scale factor of the universe. To recover the present-day Hubble constant Ho it 
is required that 

n + n R + n A = 1 , (10) 

where the terms represent the densities of matter, curvature, and a cosmological constant relative 
to the critical density of the universe. For f^A = 0, the E(z) function simplifies significantly 

E{z) = {l + z)(l + £lz) 1 ' 2 . (11) 

We assume A = 0, h = 0.75, and = 1 for the cosmological parameters of our standard model, 
motivated by the remarkable match between theory and observation of the cosmic microwave 
background and large scale structure power spectrum (e.g., Gawiser fc Silk 1998] and references 
therein). 

The SN-induced CGB depends sensitively on the evolution of the cosmic star formation rate. 
Recent progress in this field now provides measurements of this evolution to redshift well past z = 
1. We already introduced (eq. 3) the local rate density derived from H a surveys of nearby galaxies 
(e.g., pallego et al. 1995| ) to which we normalize the star formation history. Conversion of the 
H Q luminosity function of galaxies in the local universe can be related to the star formation rate 
under the assumption of an IMF. The value from pallego et al. 1995 , which we employ in this 



study, is based on the Salpeter IMF. Recent UV observations of |Treyer et al. 1997| and H a data of 
esse Sz Maddox 1998j suggest a factor two increase by redshift z ~ 0.2, which suggests a ~ 3—4 
for an evolutionary law of the form SFR(z) oc (l+z) a , i.e., a rather rapid increase of the SFR with 
redshift. 
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How much do we know about the star formation rate density as a function of redshift? 
Besides chemical evolution evidence for a significant increase in the star formation rate of the 
Milky Way as a function of look-back time, similar evidence for enhanced past star formation was 
found for faint galaxies at redshift beyond z ~ 0.3 ( Ellis et al. 1996| ; Cowie et al. 1997; Lilly et 



al. 1996). A strong increase in the comoving star formation rate density has also been predicted 



by cosmic chemical evolution models (e.g., Pei Fall 1995 ) addressing the observations of Lyman 



absorption systems in QSO spectra. The advent of the Hubble Deep Field, HDF ( [Williams et 



al. 1996 ) allowed photometric surveys to probe to z ~ 5, and Madau et al. 19*96] showed that 



the SFR appears to peak around z = 1.5, and further suggested that the SFR slowly declines to 
present-day values by redshift z ~ 5. Since then, analysis of individual galaxies in fields flanking 
the HDF ( |Guzman et al. 1998 ) has supported the rapid increase in the SFR to redshift near 



unity, and combined HST and ground-based IR photometry of the HDF (Connolly et al. 1997) 
confirmed that the SFR increases by at least an order of magnitude to z ~ 1 and that a peak in 
the comoving rate density occurs between z ~ 1 and z ~ 2. 

Although the detailed shape of the function describing the cosmic star formation history 
has not yet been determined, especially its behavior at large redshift, at least one robust and 
consistent result appears to emerge from the data: The comoving star formation rate density has 
decreased by about one order of magnitude since its peak at z p ~ 1 — 1.5 (|Guzman et al. 1998 ). 



We thus represent the normalized SFR function by a simple function (see |Yungelson &: Livio 1998 
for a similar approach) 

logMz)) = I A l0S ^ + ^ f ° r Z ~ Zp (12) 

| A log(l + z p ) — B (z — z p ) for z > z p 

where we select a "standard" case (A, B, z p ) = (4, 0.25, 1) for presentation in the figures. For 
fixed B = 0.25, we explore variations in A and z p to constrain the location of the peak of the star 
formation history with the CGB, and we also consider cases with constant star formation rate 
before the peak (B = 0) in order to address the possibility that some significant fraction of the 
star formation activity in the early universe could be hidden by dust (e.g., Hughes et al. 1998). 



In all cases star formation is assumed to be zero beyond z = 5. As far as the nuclear yields are 
concerned we apply standard supernova models, neglecting corrections due to the changing metal 
content in the host galaxies. These effects are expected to be of second order, but might be worthy 
of further study as our understanding of the CGB improves. Supernova surveys suggest a separate 
treatment of SNI and SNII, and a delay between these two event classes. We discuss this point 
below in more detail. 

The distribution of events with redshift broadens the monochromatic 7-ray lines into an 



observable continuum. We take this into account following standard prescriptions (Weinberg 1972 
Peebles 1993] ) and calculate the differential flux per unit energy. 
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3. Long-Lived Isotopes 

For isotopes that decay with a lifetime exceeding a few years the expanding atmosphere no 
longer provides enough optical depth to alter the line spectrum. The emerging lines are thus 
treated with the formalism described above, assuming that all photons have a 100% escape 
probability. We consider 26 Al, 44 Ti, 60 Co, and positrons. 

26 Al has a half-life of 7.2 x 10 5 years, and predominantly decays into an excited state of 
26 Mg via /3 + -decay and electron capture. A single 7 ray at 1.8 MeV is emitted (See Fig. 3.6 of 
Arnett 1996p . We assume that each Type II supernova ejects M26 = 1.0 x 10 4 M Q of radioactive 



aluminum ( |Timmes et al. 1995| ). 

44 Ti is believed to be the dominant nucleosynthetic progenitor of stable 44 Ca ( [Bodansky. 



Clayton, Fowler 1968| ; [Arnett 1996] ) . 44 Ti decays (half-life of ~ 60 years) to 44 Sc. The decay of 



44 Sc (half-life of 3.93 hours) to Ca generates a 7-ray photon (98.99 % of the time) at E = 1.157 
MeV. We assume a typical supernova to eject M44 = 5.0 x 10~ 5 Mq ( [Tunnies et al. 1996j ). 60 Co 
decays with a half-life of 5.3 years to 60 Fe, emitting two 7-ray photons at E = 1.17 MeV and E = 
1.33 MeV. A characteristic ejecta mass (M 60 ) is 3 x 1O~ 5 M ( [Tunnies et al. 1996j ). 



The diffuse 511 keV glow of the Galaxy (e.g. Prantzos 1993| ) is assumed to be similar to that 
other galaxies. We include this line by scaling it to the SNI rate, assuming that 3% of all 56 Co 
positrons escape from SNI and find their way into the ISM. The resulting cosmological 511 keV 
feature is shown in Figure t| Such a fraction would explain most of the Galactic annihilation line 
(e.g. Purcell et al. 19*9^ ) . However, this line will always be dwarfed by the contribution from 



annihilation of 56 Co positrons in the supernova ejecta. Averaged over the supernova event, roughly 
40% of the 511 keV photons from the 97% of the positrons that annihilate in the SNI envelopes 
escape the ejecta, which will necessarily exceed the diffuse emission corresponding to 100% of the 
photons from 3% of the positrons. 



4. Isotopes with short lifetimes 

For isotopes with short half-lives, Compton scattering in the expanding supernova is 
important and detailed 7-ray transport studies are required to determine the emerging spectrum. 
The formalism presented above only needs a small modification, through the introduction of the 
emerging differential spectrum normalized to the number of primary 56 Ni nuclei. 



4.1. Type II 

The light curve of Type II supernovae (SNII) 
from the decay chain 56 Ni — > 56 Co — > 56 Fe. Some 



supernovae 

is partially powered by the energy deposition 
of the photons emitted in this process escape 
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and thus contribute to the CGB. 56 Co ( f]/ 2 = 77.12d) decays to 56 Fe, emitting 7-rays at 0.847 
MeV (100%), 1.04 MeV (14%), 1.24 MeV (68%), 1.77 MeV (16 %), 2.03 MeV (12 %),2.6 MeV 
(17%), and 3.24 MeV (12.5%) among other lines (See Fig. 13.4 of [Arnett 1996] for a simplified 
decay scheme). Compton scatterings degrade these photon energies, leading to a broad continuum 
spectrum that develops underneath the line spectrum. A typical SNII line photon has a small 
escape probability of order 1%, because of the massive hydrogen envelope. Our treatment of 
photon transport in an expanding supernova is described in |The et al. 1990 . As standard input 



for the expanding envelope we use model W10HMM ( pinto fc Woosley 1988bj ), time integrated, to 
provide a SNII template for 7-ray photons resulting from 56 Ni decay. This source function, S(E) 
- shown in (Figure p]), is used to calculate the CGB flux from the "prompt" continuum escaping 
from SNII. The number of 7-rays per unit energy as a function of energy is thus given by 

N y (E) = N 56 S(E) (7/keV) . (13) 

The differential CGB flux is 

^I^lL{ 1 cm- 2 keV- 1 S - 1 sr- 1 ) = F Q C con , (14) 

where Fq is the same as before, but the integration over the cosmic star formation history is 
modified to 

C con = / dz E(z)- 1 V (z) S(£ 7 x (1 + z)). (15) 
Jo 

The factor (l+z) -1 in the previous C function is now absent, because of a compensating factor 
(l+z) due to the compression of the bandwidth in photon energy. 

To estimate the parameter £$N that appears in the normalization Fq we assume a Salpeter 
IMF in a mass range between 0.1 M and 125 M . If Type II supernovae result from stars with 
masses exceeding ~ 8 M , one finds £// = 0.007 (see eq. 3 in |Madau 1998 ), which is in agreement 



with arguments based on average Galactic chemical evolution ( [Timmes et al. 1997 ). For the 



ejected nickel mass in SNII we assume a value of M_ 4 = 7.5 x 10 2 ( pinto fc Woosley 1988b|) . 



4.2. Type la Supernovae 

We now estimate the contribution from SNI, which turn out to be the dominant contributors 
to the CGB in the MeV window. The method of calculation for the emerging continuum spectrum 
is identical to that used for SNII. As a standard model we use the fully mixed version of W7 
( Nomoto et al. 1984| ), integrated over 600 days (see Figure |l| ). The figure shows the reason why 
SNI contributions to the CGB dominate: the escape fraction in the MeV regime is more than an 
order of magnitude larger than that of SNII. The normalization in terms of the ejected 56 Ni mass 
is M_4 = 5 x 10 3 , i.e., 0.5 solar masses of iron is produced in an average SNI. This should not be 
considered a rigorously uniform mass, since despite the impressive uniformity of SNI light curves 



- 8 - 



it is well recognized that SNI in fact do show a significant spread in light curves, peak brightness, 
and spectral evolution. Some of this intrinsic spread is most likely due to different nickel masses. 
The observation of sub-luminous SN 1991bg ( Turrato et al. 1996| ; Mazalli et al. 1997| ) suggested a 
56 Ni mass of 0.07 M (similar to values commonly attributed to SNII), while SN 1994D apparently 
required between 0.5 and 1.0 M Q ( Vacca fc Leibundgut 1996| ), However, we are only interested in 
the mean ejecta, because the CGB data have no bearing on any individual event. 

While there is essentially no time delay between star formation and subsequent SNII because 
of the short main sequence lifetime of massive stars, a significant delay may have to be included 
for SNI. Chemical evolution studies of the Milky Way have often been used to argue for a delay 
of more than 1 Gy, based on a break in the [O/Fe] vs. [Fe/H] distribution of stellar abundances 
(see Pagel 1995] and Yoshii et al. 1996] for recent discussions). Population synthesis models 
using various SNI progenitor schemes suggest a range of possible delays from zero to several Gy 
flRuiz-Lapuente fc Canal 1998| ; |Yungelson fc Livio 1998] |Sadat et al. 1998fc Madau 1998] ) . A 



measurement of the delay time scale would go a long way towards understanding or constraining 
SNI progenitor models. Such a measurement could be achieved with deep SN surveys, measuring 



the comoving I/II ratio as a function of redshift (e.g., Ruiz-Lapuente & Canal 1998; Sadat et 



al. 1998 ). While supernova searches now routinely detect SNI at z ~ 1 ( Perlmutter et al. 1997 
[Tonry et al. 1997 ) and the comoving SNI rate at z ~ 0.4 has already been determined by the 



Supernova Cosmology Project ( Pain et al. 1997 ), a similar accomplishment for SNII will have to 
wait for the Next Generation Space Telescope, NGST flMadau 1998; ). The present-day average 
SNI fraction of the total rate is ~ 20%, and the I/II ratio is ~ 1/3 (e.g., |Cappellaro et al. 1997 ), 



but varies with redshift due to the delay between SNI and and SNII (see Fig. 5a in Yungelson 
fc Livio 1998j ). Depending on the assumed delay time between SNI and SNII we adjust such 
that the I/II ratio at z = is fixed at the observed value 1/3. This prescription fixes the SNI 
rate evolution as a function of redshift in terms of two constants (delay time, and I/II ratio), but 
is otherwise determined by the global star formation history and its associated SNII rate (fixed 
through £/). This prescription is convenient, but not necessarily realistic. However, the lack of 
rigorous theoretical or observational guidance about the evolution of the I/II ratio leaves us little 
choice. Continued efforts to determine supernova rates as a function of redshift will eventually 
lead to a better model, which would further improve estimates of the supernova contribution to 
the CGB in the MeV regime. 



5. Results 

We first present the background due to SNI in Figure § for a few choices of the Hubble 
constant (h= 0.55, 0.65, and 0.75). It is clear that supernovae can not contribute to the CGB 
above 3.5 MeV, simply because they do not produce radioactive isotopes that emit 7-ray lines 
above that energy. The presence of strong lines leads to steps in the spectrum that might be 
observable with future 7-ray instruments. If one had a dominant strong single line, the slope of 
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the cosmologically broadened spectrum would be a direct measure of the cosmic star formation 
history. In reality the line(s) are blended with the Compton continuum, and the detailed shape 
of the spectrum can not be resolved with current technology. Still, if it is true that the CGB 
in the MeV regime is ~ 100% due to emission from supernovae, then the spectrum shown in 
Figure |2| provides an independent and unique measure of the universal star formation history! For 
comparison we also display previous results derived by [The et al. 1993| , who used a simple model 
of Fe synthesis in the universe to derive the CGB. We scale their result by a factor 2/3, in order 
to take into account estimates of the fraction of iron that is contributed globally by supernovae of 
Type la ( Thielemann et al. 1991| ). 



Next we show the contributions to the CGB from Type II supernovae (Figure ||), and include 
contributions from long-lived isotopes ( 26 A1, 44 Ti, 60 Co) as well as positrons. It is clear that these 
nuclei do not contribute substantially to the MeV background. The SNII contribution is small 
compared to that from SNI. Integrating the spectra, we estimate that the II/I ratio is ~ 0.01, as 



also pointed out by The et al. 1993. 



In Figure ^| we plot the theoretical estimates together with existing measurements of the 
CGB. The standard model is apparently able to explain the bulk of the ~ 1 MeV observed flux. 
In fact, for certain choices of model parameters 7-ray overproduction occurs, so that we are in 
principle able to constrain cosmological parameters and the cosmic star formation history. From 
the comparison of the different contributions it is clear that the MeV background is dominated by 
7-rays from SNI. However, not all questions are answered. Even if the CGB is due to SNI, we 
must explain the apparent absence of a sharp drop in the observed flux above 3.5 MeV, where 
supernovae do not emit at all. Is it reasonable to expect blazars to fill in the gap in just the 
right way to explain the nearly continuous power-law behavior of the CGB? Too little is known 
about blazar spectra in this energy range ( Bloemen et al. 1995] ; Blom et al. 1995| ), The smooth 



continuation of the SMM data from 2 MeV to 6 MeV might indicate that SNI contribution to the 
flux near 1 MeV is only a small fraction. A similar problem is apparent at lower energy (few 100 
keV), where Seyfert spectra fall rapidly and supernovae create only a flat spectrum. 

Let us place this result (Figure 0) in perspective: the most plausible origin of the CGB in the 



soft 7-ray region (< 0.5 MeV) appears to be a sum of unresolved Seyfert galaxies (e.g. Zdziarski 



et al. 1993 ; Zdziarski 1996| ), which is consistent with balloon observations (e.g. Kinzer et al 



1978 ). In the hard 7-ray region (> 10 MeV), unresolved blazars are successfully invoked (e.g. 



Stecker fe Salamon 1996| ) to match high-energy EGRET observations QCheng and Mukerjee 19*98 



Sreckumar et al. 1998|) . The spectral properties of Seyfert galaxies and blazars leave a window 



around 1 MeV in which their integrated fluxes would be far below the observed level (Kappadath 



et al. 1996 ; Watanabe et al. 1997| ). It is perhaps reassuring that our standard model for the 



supernova contribution is within a factor two of the observations, which argues for SNI as the 
primary explanation of the bulk, or all, of the MeV background. This fact is exploited below 
to generate constraints on the cosmic star formation history. However, we note the discrepancy 
between the combined Seyfert/Ia prediction and the data around 400 keV, perhaps hinting at a so 
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far undetected population of sources or else suggesting that a fraction of Seyferts may have harder 
spectra than those found so far with the Compton Observatory. 

As apparent in Figure || the line features at 847, 1238, 1770, 2030 and 2599 keV due to 7-ray 
line emission from the decay chain 56 Ni — > 56 Co — ► 56 Fe are not washed out completely, but still 
imprint a step structure in the CGB spectrum. So far one (unsuccessful) attempt was made to 
detect these edges flBarthelmy et al. 1996Q , leaving a significant challenge for other current and 
also for future observations in the MeV regime. 



5.1. Progenitors of Type la Supernovae 



A satisfactory progenitor scenario for SNI has not yet emerged (e.g., [Nomoto et al. 1997 



Trimble & McFadden 1998) 


1973 




[ben & Tutukov 1984) 



1984), i.e., distinguished by the nature of the white dwarf companion. Models of the explosion are 
sensitive the mass of the white dwarf. One refers to Ch or sub-Ch models, depending on the mass 
of the unstable white dwarf relative to the Chandrasekhar limit. We might learn something about 
the right progenitor systems through observations of the delay between star formation and la 
occurrence, which one might be able to derive from the observation of the relative I/II supernova 



rate as a function of redshift (Madau 1995 ; Sadat et al. 1998; Ruiz-Lapuente & Canal 1998). 
Timescales of SNI explosion in both model classes depend on various parameters, such as the 
initial separation and the mass ratio of two WDs in DD models and length of time spent filling 
the Roche lobe in SD models ( Madau 1998| ). This time delay of SNI explosion after the global 
SFR can be ~ 10 9 yrs or as short as ~ 10 7 yrs (e.g., Ruiz-Lapuente fc Canal 1998|) , depending on 



the particular scenario under consideration. Following [Madau 1998| we employed a characteristic 
explosion (delay) time scale (17), defining the explosion probability per white dwarf, and assume 
tj to be epoch independent. Given the vast uncertainties in the la heritage, we consider three 
cases:?/ = 0.0, 1.5, and 3.0 Gy. Figure || shows the effects of these choices. Note that we keep the 
ratio of I/II = 1/3 at z=0 for any cases. 



5.2. Constraints on the Cosmic Star Formation History 

Figure ^ shows the evolution of the supernova rate as a function of look-back time (which 
depends somewhat on our choice of cosmological parameters). Figure shows the same functions, 
but displayed in redshift space. The normalization is i?//(z=0) = 1. As discussed above, the 
function for SNII is identical to (our fit of) the observed star formation rate history, while the SNI 
curve shows a delay, but is normalized to give a present-day ratio of SNI/SNII = 1/3. This ratio 
decreases with redshift (Figure ||), resembling the trend predicted by population synthesis models 
(e.g., see Figure 5 in Yungelson fc Livio 1998 ). 
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Figure |9] shows multiple curves indicating how the CGB grows with maximum redshift 
included. Here redshifts up to 0.5, 1.0, 1.5, 2.0 and 5.0 are plotted, showing that most of the CGB 
comes from SFR at redshift less than ~ 1.5, so that our results are insensitive to the possibly 
large uncertainties in measurements of the cosmic star formation rate at high redshifts. We also 
considered star formation histories in which the rate remains constant between z = z p and z = 
5, but this affects only emission from high redshift so that the correspondingly shifted supernova 
spectra contribute in an energy regime that is completely dominated by other sources. In other 
words, measurements of the MeV background are not able to constrain the SFR at high redshift. 
The situation is far better for the recent SFR history, which is the last point of discussion. 

For a given delay time r the supernova-induced CGB is mostly determined by the location 
and strength of the peak in the cosmic star formation history. In general, a shorter delay and/or 
smaller z p will increase the CGB. If the delay is zero, the CGB is inversely proportional to the 
location of z p (in our simple representation of the SFR with two power laws). Placing the time of 
the peak of star formation rate at larger redshift reduces the CGB. Of course, the total integrated 
flux also depends on the assumed I/II ratio, and the overall amplitude of the SFR. The current 
observational status of the SFR history (e.g. [Lilly et al. 1996 ; |Connolly et al. 1997 ; (Sadat et al. 



1998| ; [Madau 1998| ; |Guzman et al. 1998|) suggests values around (z p , rj(z p ) = (1—2,10—30). These 



values do not lead to a conflict with the measured CGB. Any increase in nickel yields in an average 
SNI or a global increase in the I/II ratio would eventually conflict with the CGB observations, 
which therefore provides an independent constraint on the global properties and cosmic evolution 
of Type la supernovae. While the constraints on cosmological parameters (Ho and O) do not lead 
to improvements over other methods, the CGB provides a unique handle on the unknown delay 
parameter 17. 

While the increase of the SFR with redshift by a factor 20±10 to the peak location z p remains 
a robust conclusion of recent studies of the SFR, the normalization has been challenged. A new 
determination of the local volume-averaged star formation rate from the 1.4 Ghz luminosity 
function of star forming galaxies ( Serjeant et al. 1998| ) implies a local SFR density 2— 3x larger 
than the Gallego et al. H a estimate. Tresse and Maddox (1998) have shown the H Q luminosity 
function at z = 0.2, and their data suggest a SFR twice that determined from UV measurements. 
A similar shortfall of UV-based estimates in comparison to those derived from Balmer lines was 
reported by Hlazebrook et al. 1998 who used J-band IR spectroscopy of a redshift selected sample 
of 13 galaxies from the Canada-French Redshift Survey (CFRS) to measure the H Q luminosity 
function at z = 1. These observations also indicate elevated SFR values (factor 2—3) relative to 
values derived from the UV. We already mentioned the deep submillimeter survey of the HDF 
( Hughes et al. 1998| ), which implies a star formation rate for z = 2— 4 that is five times higher than 
that derived from optical and UV observations of the HDF. The interpretation of these differences 
involves the well known fact that star formation occurs in dense molecular regions which hides the 
optical and UV emission due to substantial extinction. That much of the cosmic star formation 
activity occurs in very dusty regions is also supported by observations of the IR background in the 
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140—240 [im region, which was recently detected by the DIRBE and FIRAS experiments aboard 
COBE (Dwek et al. 1998 and references therein). The IR background data also show that the 
UV and optically determined star formation rates fall short in producing the IR background, and 
specifically require the peak star formation rate (at z ~ 1.5) to be larger by at least a factor two 
( Dwek et al. 199£| ). These various recent claims strongly suggest that the SFR(z) function we 
use as "standard model" should be multiplied by a factor 2—3 at all redshiftsl This would simply 
mean that we have to multiply the CGB fluxes by the same factor. Without the delay of SNI 
explosions this would yield a CGB spectrum in excess of the observed values, while for a 3 Gyr 
delay the model matches the observed flux. With the revised SFR values we can thus explain 
all of the CGB with emissions from Type la supernovae, and one does not need to invoke a new 
source population. However, a delay of 3 Gy is on the extreme side of the suggested values. It is 
clear that the CGB significantly constrains the properties of SNI and possible further increases in 
SFR values. If future observations can provide an accurate functional form for SFR(z), the CGB 
can be used to constrain supernova models. However, if SFR(z) estimates continue to be revised, 
the CGB provides a useful upper limit. 



6. Conclusions 

We calculated the contribution of supernovae to the cosmic 7-ray background (CGB). 
Following [The et al. 1993 we used models W7fm and W10HMM as source templates for SNI 



and SNII, respectively. Our approach differs from that of [The et al. 1993| through the use of 
the observed star formation history of the universe, obtained from Lya QSO absorption studies 
(e.g., Pei fc Fall 1995| ), galaxy redshift surveys (e.g., Lilly et al. 1996| ), broad band photometry of 



galaxies in the Hubble Deep Field (e.g., Madau et al. 1996| ), and other recent work in this very 



active area of observational cosmology. We consider time delays of 0—3 Gy between SNII and SNI 



explosions. Our estimated background spectra are similar to those derived by The et al. 1993 . We 
confirm their finding that SNII contribute very little (~ 1 %) to the CGB. SNI on the other hand, 
could explain most or all of the observed CGB spectra, for certain parameter choices. However, 
the standard model does not match the observed flux of the CGB, which suggests that either 
the current flux measurements are still an overestimate, that there could still be an unrecognized 
source population making a substantial contribution to the MeV background, or that the cosmic 
star formation rate is significantly higher than commonly assumed. This discrepancy would 
become even more serious if the SNI delays were much larger than 1 Gy, and/or if most of the 
cosmic star formation activity occurred at redshifts past z = 1. For example, combining Madau's 
rates and a 2 Gy delay leads to a CGB flux that falls short of the observations by more than a 
factor three. On the other hand, short delays (less than 1 Gy) combined with a star formation 
history that has an increase of 10-30 by z p = 1 yield a CGB flux that is just sufficient to explain 
the observations in a limited range of photon energies. We conclude from this, that the currently 
favored scenarios of SNI progenitors and their cosmic rate evolution underpredict the CGB. Recent 
upward revisions of the cosmic star formation history at all redshifts increase the predicted CGB 
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fluxes, compensating the shortfall. However, an increase by a factor 2—3 overproduces the CGB 
unless the SNI delay time scale is much larger than 1 Gyr. 

While it can not be proven, it also can not be ruled out that SNI could provide most or 
all of the CGB. If that were the case, theory predicts strong spectral steps in CGB due to line 
features from the 56 Ni — > 56 Co — > 56 Fe decay. These have not yet been observed, but remain 
a challenge to current and future 7-ray experiments. Some fraction of the CGB (mostly below 
each of the major line energies) must be non-isotropic, because the high energy end of each line 
is due to emission from the nearest galaxies, which are known to have a very non-uniform spatial 
distribution to distance of ~ few 100 Mpc (or redshift of z ~ 0.1). The angle averaged spectrum 
out to z = 0.1 integrated over energy (thus the number of photons cm~ 2 s^ 1 ster^ 1 ) is about 2% of 
the total (integrated to z = 5). Thus, to see the anisotropy in this energy range one would need 
new detectors that can detect the CGB to 1% accuracy. 

The 7-ray background in the MeV regime provides valuable constraints on global iron synthesis 
in the universe, the global rate of star formation, and yields and lifetimes of SNI progenitors. 
These constraints are complimentary to other astronomical methods. Current observations of the 
CGB in the MeV range can largely be explained as the unresolved superposition of 7-rays emitted 
by Type la supernovae. At present there is no unsurmountable conflict between the estimated 
SFR, I /II supernova rate ratio, nickel yields, la explosion time scales, and the observed CGB flux. 
However, the parameter ranges are significantly constrained and there are open questions about 
the quality of the CGB fit from the combination of Seyfert galaxies, SNe, and blazars. 
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Fig. 1. — The time integrated 7-ray continuum flux per 56 Ni nucleus for the SNII template model 
W10HMM (dashed line), and for the SNI template model W7fm (solid line). SNI dominate in the 
MeV regime, because of their larger fraction of escaping 7-rays. SNII dominate at lower energies, 
but in this regime the CGB is overwhelmingly determined by X-ray emission from Seyfert galaxies. 
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Fig. 2. — The CGB contribution from SNI based on the SFR evolution function r\ (eq. 12) with A 
= 4, B = 0.25, z p = 1, and 1.5 Gy time delay for SNI. Three values for the Hubble constant are 
presented. For comparison we also show the results of The et al. 1993| , multiplied by 2/3 (based on 
the assumption that 2/3 of the cosmic iron production can be attributed to Type la supernovae; 
Thielemann et al. 199l| ). 
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Fig. 3. — Type II supernovae contributions to the CGB from 56 Ni — > 56 Co — > 56 Fe decay 7-lines 
and their corresponding comptonized 7-ray photons including long-lived line afterglows due to 26 Al, 
44 Ti, 60 Co, and positrons. 
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+ - HEAO-1 A4 MED (Kinzeretal. 1997) 
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Fig. 4. — The CGB from SNI, SNII (including long-lived isotopes), in comparison to recent CGB 
measurements from HEAO-1 A4 flKinzer et al. 1997|) , COMPTEL ( [Kappadath et al. 1996|) , and 
SMM ( Watanabe et al. 1997| ) . For the SMM data the dashed-dotted lines indicate the ± la region 
of uncertainty of the averaged spectrum. Predictions for the background contributions from Seyfert 
galaxies (Zdziarski et al. 1993) are also shown (dashed line). 
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Fig. 5. — Effects of the characteristic explosion/delay time scale (tj, see text) of SNI relative to 
SNII which directly trace the cosmic star formation history. Three cases are shown:r/ = 0.0, 1.5, 
and 3.0 Gy. Note that we keep the SNI/SNII ratio at z = fixed at 1/3 in all cases. The comparison 
to the preliminary SMM estimates of the CGB ( Watanabe et al. 1997|) shows that large delays lead 
to a CGB flux that significantly falls short of the observed spectrum. 



-23- 




Fig. 6. — Normalized SN rates in function of look-back time. The SNII rate (dashed line) is 
normalized at today. The SNI rate (solid line) has a delay time from a birth of a WD to the 
SN explosion. For simplicity we adopted a mean lifetime of SNI progenitors of 1.5 Gy which was 
obtained from the chemical evolution study ( |Yoshii et al. 1996 ). Note that different £sn i n Eq. (6) 
are used for each type of SN. 
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Fig. 7. — Same as Figure g but function of redshift (z). Transformation between z and t 
involves a model dependent look-back time. 




Fig. 8. — The ratio of the SNI rate to the SNII rate as a function of redshift. The ratio is nearly 
constant for z > 1 but rapidly increases towards small z, reaching 1/3 at z = according to our 
adopted normalization for the present-day derived from local supernovae surveys (see text). 



-26- 




Fig. 9. — Curves indicating how the CGB grows with integration to maximum redshifts 0.5, 1.0, 
1.5, 2.0 and 5.0 (no SNI delay was assumed for this calculation). The figure illustrates that most 
of the CGB is due to emission from redshifts smaller than the peak resdhift of the cosmic star 
formation history. 



